The spectral properties of noncondensate particles in Bose-condensed atomic 

hydrogen 
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The strong spin-dipole relaxation, accompanying BEC in a gas of atomic hydrogen, determines the 
formation of a quasistationary state with a flux of particles in the energy space to the condensate. 
This state is characterized by a significant enhancement of the low-energy distribution of non- 
condensate particles, resulting in the growth of their spatial density in the trap. This growth leads 
to the anomalous reconstruction of the optical spectral properties of non-condensate particles. 



The discovery of Bose-Einstein condensation (BEC) in 
a trapped dilute gas of atomic hydrogen has opened 
a new page in the studies of metastable gaseous systems 
at ultracold temperatures. For a long period, trapped 
atomic hydrogen remained a very promising system to 
achieve BEC (see, e.g., However, only after the in- 

volvement of the forced evaporative cooling mechanism 
typical for the works resulted in the BEC discovery in al- 
kali metal vapors (^-Q] the achievement of BEC in atomic 
hydrogen became a reality. 

The experiments reported in have brought out a 
set of specific features inherent in the hydrogen system. 
The anomalously low rate of a three-particle recombi- 
nation allows one to reach the record density of a Bose 
condensate of n c ~ 10 16 cm~ 3 with the total number of 
atoms N by several orders of the magnitude greater than 
that in the experiments on alkali gases. The lightest mass 
of hydrogen atoms in combination with so high density 
determines the superfluid phase transition temperature 
T c exceeding that for the experiments with alkali gases 
by two orders of the magnitude. 

One of most interesting features of the hydrogen exper- 
iments is impossibility to achieve a relatively high concen- 
tration of the condensate fraction. This result is a direct 
evidence for the phenomenon of the "burning of a conden- 
sate" H related to a strong increase of the gas density at 
the condensation and, as a result, drastic enhancement of 
the spin relaxation rate (see also f|). The low rate of the 
evaporation cooling due to an extremely small scattering 

o 

length a w 0.65A makes this phenomenon responsible for 
the kinetics of the formation of a condensate in a gas of 
atomic hydrogen. This led to the prediction ||| that the 
concentration of the condensed fraction cannot exceed 
several percents. It is worthwhile to note that the spatial 
density of the condensate may be large. 

In the experiments E| an unexpected phenomenon has 
been observed. The shift of the line of the ls-2s tran- 
sition for non-condensate particles is found to increase 
abruptly after the BEC transition. This is especially sur- 
prising since the condensate volume is only about 1% of 
the volume of the thermal cloud and the estimations pre- 
sented in (!]] could not explain the observed picture. The 
aim of the present paper is to reveal the nature of this 
interesting phenomenon. 



Under conditions of the experiments the interpar- 
ticle collision time r is small compared with the lifetime 
of the system rr,. As a result, a quasistationary state 
sets in for r < t < tl . For T > T c , this state is close 
to the equilibrium one. However, for T < T c the intense 
losses of the particles from the condensate due to spin re- 
laxation determine the appearance of the quasistationary 
state with the compensating flux of particles into the re- 
gion of low energies. As turns out, the distribution func- 
tion for the non-condensate fraction of the gas in such 
state differs drastically from the equilibrium distribution 
function. An important feature of this "steady-flux" dis- 
tribution is a sharp increase of the number of particles 
in the low-energy range of the spectrum. For the trap 
geometry, this results in a strong increase of the den- 
sity of normal fraction and therefore causes the anoma- 
lous transition line shift for non-condensed atoms. Note 
that the attempt to explain the anomalous increase of the 
non-condensate density by introducing a rather artificial 
model of the formation of metastable dense droplets has 
recently been undertaken in [ ^0[ . 

We confine ourselves to the time scale of r <C t <C tl 
and suppose that the quasistationary regime sets in. In 
order to simplify the analysis we will consider a gas in 
a spherically symmetric harmonical potential with fre- 
quency u>. We suppose that the condensate is in the 
so-called Thomas-Fermi regime at T < T C) that is, 
n c (0)Uo 3> huj where Uq = Airfi 2 a/m and m is the atom 
mass. This means that the condensate density is deter- 
mined by the expression 



n c (r) = 



2/i c — muj 2 r 2 

2Un 



(1) 



valid for r < r M where r M = (2/Zc/mw 2 ) 1 / 2 and (i c = 
n c (0)£/g is the chemical potential of the condensate. 

At T < T c the loss of the particles is determined mostly 
by the condensate 
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h c {r) = --L^n 2 c (r). 



(2) 



Here is the spin relaxation coefficient determined 
for a normal gas and the factor of 1/2 appears from the 
two-particle density correlator for a condensate. Note 
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that, in spite of the spatially inhomogeneous loss of Bose- 
condensed particles, the condensate retains Thomas- 
Fermi density profile. The point is that the redistri- 
bution of particles in the condensate, induced by the 
mean-field interaction, is much faster than their loss. In- 
deed, the characteristic evolution time in the condensate 
is t = h/n c (r)Uo, while the characteristic time of the 
particle loss is t; = 1/ L^n c {Q). In the atomic hydrogen 
ti/t ~ 10 5 . 

With the condensate density distribution (|l|) we find 
for the total particle loss rate 



Q = 

* 105 



(3) 



Let us consider the kinetic equation for the normal frac- 
tion in a harmonic potential. We are interested in par- 
ticles with energies e > a 3> tku. In this case we can 
disregard the discrete structure of the trap energy levels. 
Under the assumption of ergodicity the quantum Boltz- 
mann equation acquires the form 



n(ei) = I co ii{e\, [n]). 



(4) 



where n(e) is the statistical average of occupation num- 
bers. Within the energy interval where the occupation 
numbers are large 
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This kinetic equation differs from that for a homogeneous 
case (see, e.g., p"l| , |l2|]) only by the exponent in expres- 
sion (H) (see, e.g., ]13| , p^ |), which is characteristic for a 
harmonic potential. As we will see below, the increase of 
spatial density, accompanying BEC, comes from the en- 
ergy region where n(e) 3> 1. We took the latter condition 
into account explicitly, deriving Eq. (||). 

The density in the energy space can be written as p e = 
g(e)n(e) where g(e) = e 2 /2(hid) 3 is the density of states 
in the harmonic potential. We can rewrite Eq. (Q) in the 
form of the continuity equation for the energy space 

dt de K ' 

with the flux Pie) 



P(e) = f de'g(e')I coll (e', 
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(8) 



In the quasistationary regime P(e) = Q = const. This 
means that the flux <JsJ> has the same value for an arbi- 
trary e. Taking into account Eq. (||), this requirement 
can be fulfilled if the distribution function has the form 



n(e) 



A 



(9) 



where £ = e/huj. Let us substitute this distribution func- 
tion into Eq. (||). Introducing dimensionless variables 
= £j/ei, after simple transformations we find 



where 

ioh) - 



I coll = TU(M 2 A 3 / ( 7 )£l- 37 



(a+c 3 7 -c 2 7 -i)(6^4)- 7 



Inserting these relations into Eq. 
flux 



|), we find for the 



P(e) oc 



4(7) 
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The requirement of the independence of the flux on e 
leads to 



7 =3 



(10) 



and simultaneously to ^0(7) = 0. The last result is fairly 
evident: I co u is equal to zero in a stationary case. At the 
same time the derivative I (j) is finite and provides the 
permanent flux P = Q to the low energies region. The 
obtained results are in correspondence with the general 
analysis of Zakharov (see |ll| and also ) , showing that 
the equation I co u = has a nontrivial solution relevant 
to a steady-state particle flux in the energy space. 
For the relation between A and Q, we obtain 



A = C 



( Q 



\w(huj) 



1/3 



(ii) 



One can estimate the dimensionless numerical coefficient 
C using the total number of non-condensate particles N nc 
if one assumes approximately that Eqs. (|9|),(^o|) remain 
valid up to the maximum energy. 

It is worth to note that the numerical calculations of 
the BEC kinetics demonstrate the formation of the dis- 
tribution of the form (with 7 > 1) before the real 
growth of condensate starts (see, e. g., [|l2 15 1). 

Using n(e) from Eq. ([)]), one can find the spatial den- 
sity distribution n{r) for the normal fraction. The direct 
numerical calculation shows that the main contribution 
to n(r) comes from the energy levels fi <C £ < T. These 
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levels are weakly affected by the presence of the conden- 
sate (the volume occupied by the wavefunctions of a par- 
ticle is much larger than the condensate volume) and by 
the interparticle interaction in the thermal cloud. This 
means that within a reasonable approximation we can 
neglect the renormalization of the levels and use oscilla- 
tor vavefunctions "0p( r ) where p = (n r , I, m). The energy 
region (e — /i) ~ fj,, where the modification of the level 
structure is noticeable, gives a negligible contribution to 
the spacial density, and we ignore this modification. Un- 
der the assumption of ergodicity 



n(r) =^|V P (r)| 2 n(e E 



(12) 



Since the spacial density of the thermal cloud is weakly 
sensitive to the low energies, at the presence of the con- 
densate we simply truncate the sum (|l2| ) by the condition 
e p > Mc- 

Since the concentration of the condensate fraction is 
confined to several percents, the main part of the par- 
ticles at T < T c is in the non-condensate fraction and 
N nc N. At T > T c the distribution function for low 
energies at the actual absence of the flux compensating 
losses has an equilibrium form n(e) = T/(e + \fx\). Com- 
paring this distribution with Eqs. (PJlOl), one can see that 
the maximum density of the non-condensate fraction, de- 
termined in the trap mostly by the low-energy region in 
the presence of flux at T < T c , can exceed significantly 
the equilibrium value of the density at T > T c . Note that 
in the presence of condensate the thermal equilibrium is 
absent and therefore the temperature may be regarded 
only as a characteristic energy scale for the system. 

For the Doppler-free two-photon excitation, the spec- 
tral distribution is determined by the red shift of the 
absorption line. The shift is caused by the change of the 
coupling (scattering length) of an excited atom compared 
with the atom in the ground state. In the quasiclassi- 
cal approximation the shift is proportional to the local 
density of particles, coinciding with n(r) beyond the con- 
densate region. Thus, the appearance of a condensate, 
accompanied by a sharp growth of n(r) for small r, causes 
a significant increase of the shift of the non-condensate 
spectral distribution. 

For an optically thin sample, in the local density ap- 
proximation the Doppler-free spectral distribution for the 
two-photon absorption is proportional to the density dis- 
tribution G(n). For the spherically symmetric configura- 
tion, this distribution reads 



G(n) = Anr 2 (n) 



dr{r 



dn 



(13) 



For 



where r(n) is determined from condition n(r) = 
the non-condensate fraction, n(r) is determined by (12) 
Below we present approximate quantitative estima- 
tions for the described picture on the basis of the ex- 



perimental data 0,0]. Since in the experiments the col- 
lision time r ~ 10~ 3 s is much less than the lifetime of 
the system tl ~ 5s, the quasistationary regime sets in 
for times r < f C tj, as explained above. The maximum 
condensate density is found to be n c (0) ~ 5 x 10 15 cm~ 3 . 
We will consider a spherically symmetric harmonic trap 
with oj = (cj^cjj) 1 / 3 = 3.4 x 10 3 Hz. Here uj_ and lu z 
are the frequencies of the cylindrical trap used in 

For L< 2 ) = 1.1 x 10 15 cm 3 /s and a = 0.65A [UH], we 
find Q = 1.9 X 10 9 s _1 for the particle loss rate. In the 
Thomas- Fermi approximation, the number of particles in 
the condensate for the given parameters is N c ss 1.2 x 10 9 . 
The largest uncertainty in the experimental data jlj is re- 
lated to the value of the relative condensate fraction \- 
Here we use some average value of x ~ 0.07. Then the 
total number of particles is N « 1.7 x 10 10 . The critical 
temperature of the BEC transition is T c w 60/iK in this 
case. 

Let us assume that non-condensate particles are con- 
centrated within the energy interval [/*,£*], having dis- 
tribution (H). Comparing the energies of a system with 
the distribution (||) and with the equilibrium distribu- 
tion at temperature T, we find e* = 4.7T C (T/T C ) 4 for the 
same number of particles in the both cases. The param- 
eter e* ^ /i. The increase of the spatial density of non- 
condensed particles originates from the occupation of the 
energy interval where e is significantly smaller than e*. 
The results are only weakly sensitive to the behavior of 
n(e) near the upper bound of the energy spectrum. Our 
approximate definition of e* is made with the aim to com- 
pare the results below and above T c , assuming the same 
number of particles. Taking into account the estimated 
values of N nc , Q and W from Eq.(|7) we can determine 
the coefficients of A and C in (0,y). At T » 0.9T C we 
obtain A « 2.9 x 10 5 and C » 2.5. 

Using distribution function (^J) with the values found 
for A and e*, one can determine the spatial density dis- 
tribution (|l2| ) for the non-condensate fraction. 
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FIG. 1. The spatial density 
non-condensate fraction (N — 1.7 x 10 lu and T c = 60/iK): 
Ti = 0.9T C (curve 1), T 2 = 70/iK (curve 2) and T 3 = 120mK 
(curve 3). 
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FIG. 2. The spectral distribution of two-photon absorption 
for the non-condensate fraction. The parameters are the same 
as in FIG.l 

The dependence n(r) is shown in FIG.l (curve 1). For 
comparison, we present the density profiles for two tem- 
peratures above T c : (curve 2) and 120/xK (curve 3). 
The length in the figure is given in units of the character- 
istic trap size Iq — y/fi/muo. The absence of condensate 
at T > T c allows one to neglect the spin relaxation losses 
and therefore to use the equilibrium distribution with the 
chemical potential fixing the same number of particles. 

Knowing n(r), we can find G(n) @ and therefore 
the spectral distribution of two-photon absorption. To 
make a comparison with the experimental results jl| , we 
use the relation Av(n) = r\n for the shift of absorption 
line and, according to Q], assume that the coefficient r\ 
has the same value both for the condensate and for the 
non-condensate fraction. Regardless of a specific value 
of the scattering length for the interaction between an 
excited atom and atoms in the ground state, one can 
use the relation Av(n) = Av c (n/n c {Q)). Here Av c is 
the experimental value of the shift corresponding to the 
maximum condensate density n c (0). 

The dependence I{Av) for the density distributions 
presented in FIG.l is shown in FIG. 2, Av being a shift 
of the laser frequency. Calculating I(Av), we took into 
account the experimental linewidth of 5v ~ 3kHz (see 

§)• 

The curves plotted in FIG. 2 reproduce the character- 
istic behavior and scale of the reconstruction observed in 
jjj for the two-photon absorption spectrum of the normal 
fraction with the appearance of a condensate. The nature 
of this reconstruction lies in the formation of the special 
energy distribution of particles due to the presence of the 



quasistationary flux of particles into the small condensate 
region compensating the spin relaxation losses. 

Note in conclusion that all the presented calculations 
have an estimational character due to adopted simplifi- 
cations and uncertainty in a self-consistent set of exper- 
imental data. 
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